Pionic Decays of D SJ (2317), D si (2460) and S sj (5718), S sj (5765) 
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We estimate pionic decay widths of the narrow charm-strange resonances D 3 j (2317) and D s j (2460) 
using the 3 Po model. Their one-pion decays occur through -n-ir mixing while the two-pion decays 
of D S j , (2460) occur through the virtual /o(980) meson. The mixing between 3 Pi and 1 P\ states 
enhances the single pion decay width of D s:) (2460) and suppresses its double pion decay width 
significantly. The two-pion decay width of D SJ (2460) is much smaller than its one-pion decay 
width. As a byproduct, we also calculate pionic decay widths of P S j(5765), B s j(5765) mesons in 
the (0 + , 1+) heavy doublet. 
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I. INTRODUCTION 

In April 2003, Babar Collaboration reported a new charm-strange state D s;) (2317) in the D s 7r° channel. Its spin- 
parity is J p = + and its mass is below the DK threshold 2.36 GeV. Later CLEO Collaboration reported another new 
charm-strange resonance £> si (2460) in the D*ir° channel with J p = 1+ below the D*K threshold Q. It's tempting to 
classify these two states as the (0 + , 1+) P-wave cs doublet. But their observed masses are more than one hundred MeV 
lower than quark model predictions |3j. There have been heated debates about its und erly i ng structure and orig in of 
its low mass in the literature [i II IE 11 11 II El El El El El El E^^ 

Although several non-conventional schemes such as DK molecules yj , four quark states [a 0J E3, Ell fl2L fl3l [lj, [l]| 

and Dtt atom mi were proposed, evidence is gradually accumulating that these two narrow resonances are ordinary 

cs states H H 03 El El HI m il m m I23. 

Another interesting property of D S j(2317) and Z? S j(2460) mesons is their extremely narrow widths. The decay 
channels £> SJ (2317) — » D K and D S j(2460) — + D* K are forbidden by kinematics. Therefore, their possible strong 
decay modes are one-pion and two-pion decays. The two-pion decay occurs via a virtual meson such as /o(980). The 
one-pion decay mode breaks the isospin symmetry and happens through r)-ir° mixing |30| : D S j(23l7) — > D s i] — > D s ir°, 



£) sj (2460) D*rj -> Z?*7r u . The 77 



mixing is described by the isospin violating piece in the chiral lagrangian 



r — 



4(m u + m d ) 



(1) 



where £ = exp(z7r// 7r ), tt the light meson octet and m q is the light quark mass matrix. Such a mixing is suppressed 
by the factor — r " d ~^|^ nd . Numerically the isospin violating effect is 0(1O~ 2 ) in the amplitude. While the isospin 

conserving strong decay width is O(10 2 ) MeV, one would naturally expect the one-pion decay width of D S j(2317) 
and £> S j(2460) to be around several tens keV. 

Although quantumn chromodynamics (QCD) is widely accepted to be the correct theory of strong interaction, our 
present understanding of the strong decay mechanism is rather limited. First-principle calculation of decay matrix 
elements on the lattice still has a long way to go. For some limited special cases, complicated tools such as light-cone 
QCD sum rules may be used to calculate the coupling constant |2q. In order to understand the vast amount of strong 
decay data, we have to turn to phenomenological strong decay models. Among them, the 3 Po model (or quark pair 
creation model) is the simplest and most successful one. In this work, we employ this model to calculate the decay 
amplitude of D S j (2317) and D S j(2460) mesons assuming they are cs states. A short review of 3 Po model is given in 
Section [n] The analysis of the one-pion decays is presented in Section ITTT1 The two-pion decay mode of D S j (2460) is 
given in Section llVI The last section is a short summary. 
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II. THE 3 P MODEL 
A. The model 

The 3 Pq model was first introduced by Micu in 1969 [3lJ and further developed by the Orsay group in the 1970s 
|32| . According to this model, a qq pair with J PC = ++ is created from the vacuum when a hadron decays. This 
created qq pair carries the quantum number of the vacuum. The new qq pair, together with the qq within the parent 
meson regroups into the outgoing mesons via quark rear rang ement process, which is shown pictorially in Fig. ^ The 
transition operator in the nonrelativistic limit reads |32l l33| 



r = -3 7 ^(lm;l - m|0 0) (2tt)* f d 3 k 3 d 3 k 4 5 3 (k 3 + k 4 ) ^i m (k 3 - k 4 ) xl%n <j)f w 34 bj(k 3 ) 4(k 4 ) (2) 

where 7 is the dimensionless pair creation parameter in the model, and indicates the strength of coupling. y" l (k) = 
k l Y l m (9k, <j>k) = — £m • k is a solid spherical harmonic function that describes the momentum distribution of the 
created pair. The meson state is defined as [35|: 

\A(n A 2SA+1 L AjAMjA )(p A )) = ^2~E~ A Y, (LaM La S a M Sa \J a M Ja ) 

Ml a Ms a 

X J d 3 k 1 d 3 k 2 5 3 (k 1 +k 2 -p A )^ nALAMLA (k 1 ,k 2 )x 1 s 2 A M SA <l>A^A\ ?l(kl)? 2 (k 2 )) 

(A(p A )\A(p' A )) = 2E A (2nfS 3 (p A - p' A ) . (3) 

The subscripts 1 and 2 refer to the quark and antiquark within the parent meson A respectively, ki and k 2 are 
the momentum of the quark and antiquark, k^ is their relative momentum and is the momentum of meson A. 
Sa = s qi + s q 2 is the total spin. J A = L A + S A is the total angular momentum. 
The S-matrix reads 

(f\S\i) = I + i(2ir) 4 S 4 ( Pf - Pl )M Mj A M ''B M ->c . (4) 
The decay helicity amplitude of the process A — > B + C in the meson A center of mass frame is 

M Mj a M Jb Mj c ( A _, B Q 

= y/8E A E B E c jY (LaM La SaMs a \JaMj a )(LbM Lb SbMs b \JbMj b )(LcM Lc ScM Sc \JcMj c } 

M La Ms a , 
M Lb ,Ms b , 
M Lc ,Ms c ,m 

x<lm;l -m| 0) (xf oMsG Xf B M 3B \xf A M SA Xf-J{^ l i\^f) C^ c (p) > 

(5) 
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where the spatial integral I m (A, BC) is defined as 



I 



M LA ,m 

M Tl „,M T . 



(p) = (27r)5 yd J k 1 d d k 2 d d k 3 d d k4^(k3+k4) ( 5 J (ki+k 2 -pA)^(ki+k3-pB)^(k 2 +k4-pc) 

x *Pn B L B M LB ( k l, k 3) C C L C A/ Lc ( k 2,k 4 ) ^„ A L A M L J)^lM) D^Oa ~ k 4 ) (6) 



We employ the harmonic-oscillator wave functions for S-wave states |3J] : 



^(Id.ka) 



i? 2 \3/4 



exp 



1 



(k x - k 2 ) 2 ii| 



where the parameter i? is the meson radius. For the P-wave state, we use [34 



(7) 



^ =i (ki,k 2 ) = zy-_y lm (k 1 -k 2 )cxp[--i?^(k 1 -k 2 ) 
The spin matrix element can be written in terms of Wigner's 9j symbol |32| 

Ms„ XS B M S „ aM s . Xl —ml 



= J2 (ScM Sc ;S b M Sb \S M.)(S M s \S a M Sa ;\ - m) [s(2S B + 1)(2S C + 1)(2S A + 1 

Sbc-Ms 




The relevant flavor matrix element for D S j(23l7) — > D s r] is 

1 



1 



(<Pc \<Pa Vo > = \~^{ uu + dd + s 3 S4}{cis 2 } {cis 3 } — {uu + dd~ 2s 4 s 2 } / ( ^ 
With the Jacob- Wick formula the helicity amplitude can be converted into the partial wave amplitude [36 



M JL {A -► BC) 



V2L +1 
2J A 



- (L0JM Ja \J a M Ja ) (J B M jB J c M Jc | JMj A )M Mj * Mj b M -'c (k fl ) 



Mj„,M Jr 



(8) 



(9) 



(10) 



where Mj A — Mj B + Mj c , J = Jb + 3c- The decay width in terms of partial wave amplitudes using the relativistic 
phase space is: 



JL\2 



8irM 



JL 



For the three body process A — > BCD, the width is 

1 1 



r 



(2tt) 5 16M^ 



|p B | \ Pc \ drnco dQ B d£l* c 



(11) 



B. Input Parameters 

Our convention follows Ref. [3^ where the value of 7 = 6.9 is \/96t: times larger than that used by others groups 
[37I l38l|. The masses of 77, /o(980) and charm-strange mesons are taken from PDG j^. Since the (0 + , 1 + ) beauty- 
strange mesons have not been discovered yet, we use the theoretical estimate of their masses from Ref. 4] in our 
calculation. For the (0 + , 1 + ) D S j and B S j mesons, we use the same parameter R from Ref. . We collect all these 
values in Table [I] 



III. SINGLE PION DECAYS OF D sj (2317) AND D sj (2460) 

Z?sj(2317) and D S j(2460) lie below the DK and D*K thresholds. They can only have pionic decay modes. The 
isospin-violating single pion decay occurs through r) — 7r° mixing as shown in Fig. [21 



4 



TABLE I: Value of mass and R in our calculation . 

Ds(Q-) £>;(!-) J> s j(2317) D^j 1 (2460) £>^(2460) B s (0~) g^Tj B a (0+) 5 s lpi (l + ) ^ Pl (l + ) V /o(980) 

mass (GeV) 1.968 2.112 2.317 2.460 2.460 5.370 5.417 5.718 5.765 5.765 0.548 0.980 

R (GeV -1 ) 1.37 1.59 1.89 1.85 1.89 1.37 1.59 1.89 1.85 1.89 1.47 2.7 




FIG. 2: Pionic decay of D S j (2317) via r) — ir° mixing. 



The decay amplitude of D S j(2317) 

M jl (d sj {2317)^ D s + tt°) = V JL (D SJ (2317)^ D s +ii) ■ — • V( V - tt°) . (12) 



The the V jl [d sj (2H7) -► D s + rj\ term is can be calculated using 3 Pq model. The r\ — tt° mixing is determined 
by the up down quark mass difference [3fH |: 

■Vin — TT = i ; —i — o Tr o r 



ml -m\ w ; 4 m s - m "+ m " 4 w " 

which vanishes in the limit of exact isospin symmetry. The isospin violating factor S^o^ — j3fj- 

In the heavy quark limit, there is one 1 + state with ji = L + s q = i + in the (0 + , 1 + ) doublet. We denote it as 
\l + ,ji = \ + )- The other 1+ state = | + ) belongs to (l+,2+) with ji = | + . ji is a good quantum number 

when rriQ — > oo. Therefore ji = L + s q = * and ji = L + s q = ^ states don't mix with each other. They are exact 
mass eigenstates. Moreover, it was shown that the mixing of these two states is small even if one considers the l/m c 
correction in Ref. 4Q]. In other words, the physical state Z? S j(2460) is quite close to |1 + , ji = \ ). 

In the quark model, the basis states are 1 P\ and 3 Pi states. Generally speaking, D S j(2460) should be the linear 
combinations of 1 P\ and 3 Pi states ED 



,3l = \ ) 
3 + 



= cos<9 | x Pi) +sin0 | 3 Pi) , (13) 



\l + ,jl = - ) = -sin6\ 1 P 1 )+cos6\ 3 P 1 ) . (14) 

The mixing angle Orq = — tan -1 y/2 in heavy quark limit |37^ . With finite charm quark mass, the mixing angle is 
not precisely known. However, 9 is not expected to differ dramatically from — tan -1 \J~2~ according to Ref. poT ]. 
The the partial wave decay amplitude and width of D s j (2317, 2460) are listed in Table ITU where 



W = (2tt)* tt-4 7<y w o„ (r 5 a R%R S c) 2 (1) 2 ^SMaEbE* exp ( - ip 2 P 



2 Rb + -^C 



yj v „-^_„ —r y jr --a y y , (15) 

V = R 2 A + R B + R 2 C . (16) 

The dependence of the single pion decay width of D S j(2317) on the mixing angle is shown by Fig. |SJ The mixing 
between 3 Pi and 1 P\ states enhance Z? S j(2460)'s single pion S-wave decay width, and decreases the D-wave width. 
We ignore the D-wave width because it is much smaller than the S-wave decay width. 
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TABLE II: The single-pion decay width (in keV) of L> SJ (2317, 2460) and B 3J (5718, 5765). 



Decay mode 
A^[Btv°] l 



Amplitude 



T(keV) 



D SJ (2317) [D 3 ir°] s 
B s (0+) -> [B S (Q- 



-W 



1 - 



- (£>> u )s 
f-B.fl 



32 
35 



L> sj (2460) 



48V3 



-12 



(— cos# + -\/2sin6 



33(fe Q ) 
38(g gQ ) 



35 



30 - 

> 

0) 25 - 



15 




FIG. 3: Variation of B s (5765)'s (dashed curve) and D s j(2460)'s (solid curve) single pion S-wave decay width with the mixing 
angle. 



IV. DOUBLE PION DECAYS OF (2460) 



Parity and angular momentum conservation forbid the double-pion decay mode of Z? SJ (2317) mesons. In contrast, 
D S j(2A60) can decay into D s inr and D*tttt final states. There is no reliable way to calculate the multiple pion decay 
width. One has to rely heavily on specific models to estimate the double-pion decay width. For example, the double 
pion decay of _D S j(2460) meson was assumed to occur via a virtual a meson in Ref. Q. 



TABLE III: The double-pion decay width (in keV) of D sj (2317, 2460) and B si (5718, 5765). The subscript A refers to D sj (2460) 
and B s (l + ) respectively, B refers to D s , Dg, B s (0~) and B s (l~) respectively, and X refers to the /o(980). 



Decay mode 
A — > [B/o(980)]l -> [Btx-k] l 




v(A-[B/o(980)]i,) 






hi 


r 2 


D 3J (2460) -> [£> s /o(980)]p -> [D s nn] P 
B s (0+) -> [B s / (980)]p -> [B s (0-)tvtt]p 




' 4 (2-5^) + (l-^)ift k B 2 ~ 


cos6> + VF'sin6> 


2.6 
0.45 


6.9 
1.0 


D SJ (2460) -> [D s */o(980)]p -» [D* s nn] P 
B s (l+) -> [B s /o(980)]p -> [B s (l-)7r7r]p 


VF'cosf + ^W" 


4 (l- 5 4)+(l-f 




sin# 


0.055 
0.044 


0.13 
0.12 



In our opinion, /q(980) is well established and less controversial than the sigma meson. Moreover, its width is 
measured to be 40-100 MeV j3s|. /o(980) mainly decays into pions. This process can be described by the effective 
Lagrangian 

£/ (980)7r7r = 9 '/ (980)inr /o[27T + ~ TT~ + 7T 7r ] . (17) 

The coupling constant gf Tnr is 0.83 ~ 1.3 GeV. In the following we assume the double pion decays occur through the 
decay chain with the help of a virtual / (980) meson 46J: D„-(2460) -> D s + / (980) -> D s + 2n and £> sj (2460) -> 
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D,j(2460) / (980) 



r(keV) 



FIG. 4: The double-pion decay modes via a virtual /o(980). 
(a) 



(b) 




r(keV) 




FIG. 5: Variation of the double-pion decay width for the process (a) D 3J (2460) — > D s tttv and (b) _D SJ (2460) — > Dlnir with 
gf Q7T7r = [0.59, 0.93] GeV and mixing angle 9. The shadow is the possible value of the decay widths. 

D* + /o(980) — *■ D* s + 2tt. We further assume that the decay orbital angular momentum L exists between /o(980) and 
D s (or D*) only. We want to remind the readers that the above working assumption is only used to make a rough 
estimate of D S j (2460) 's double-pion width. 



We first consider the decay chain D S j(2460) — > D s + /o(980) — > D s + 2tt. The decay amplitude is 
M(D sj (246Q) -^D s + 2tt) = v(p s] {2Am) ^D s + / ) ■ , 2 - _ 2 ■ y/\^ g M 



k 2 f 

JO 



(18) 



./() 



with Att+tt- = 2, A^o^o = 1. The vertex v(z^ (2460) -» D a + / ) can be obtained in 3 Po model, which is the function 

of momenta of D' s s momentum ke- The process -0,^ (2460) — > D* + /o(980) is Similar. We list the amplitudes and 
widths of double pionic decays of D S j(2A60) and i? s (5765) in Table ITTll where 



2 ~ 

(v) 2 ^ 8EaEbEx l kB l ex P 



— -kR 2 i?^ 



m 



V 



The dependence of their double pion decay widths with the mixing angle is presented in Fig. [S] 

For the pure academic purpose, we also estimate the doubly isospin- violating P-wave decay: D SJ (2460) — > 
£) S j(2317) + i] -> 1^(1869) + r] + r] -> D s (1869) + tt° + ir°. Numerically the width of the above decay mode is 
10 -6 eV, which is tiny and completely buried by the isospin conserving double-pion decay mode. 



V. DISCUSSION 



Assuming both D S j(23l7) and _D S j(2460) are ordinary cs mesons, we have calculated their strong decay width using 
the 3 P model. The single pion decays are isospin-violating and occur through rj-ir mixing. We have collected the 
single pion decay widths available in literature in Table HVI Our results agree very well with a recent light-cone QCD 
sum rule calculation 1251. 
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TABLE IV: Single-pion decay widths (in keV) of -D s j (2317) and D 3 j(2460) mesons from various theoretical approaches. 







this work 


f25] 


f42] 


m 


[43] 


[44] 


[8J [45] 


D* sJ (2317) 




32 


34-44 


7± 1 


21.5 


- 10 


16 


10-100 150 ±70 


A,j(2460) 


- D>° 


35 


35-51 


7± 1 


21.5 


~ 10 


32 


150 ± 70 



The double pion decays of D S j (2460) are allowed by isospin symmetry. But they are suppressed by three-body 
phase space. Under a rather crude assumption that such decays occurs with the help of a virtual /o(980) meson, 
we have estimated the double pion decay width to be around 2.6 ~ 6.9 keV for D SJ (2460) — * D s + 2ir mode and 
0.055 ~ 0.13 keV for D S j(2460) — >• D* + 2tt mode, depending on the total pionic width of /o(980). The double pion 
decay width of D S j(2460) is numerically much smaller than its single pion width because of the cancellation from the 
mixing of 3 Pi and 1 Pi states. 

Putting the single and double pion decay modes together, the strong decay width of D S j (2317, 2460) is less than 50 
keV. Combining the radiative decay width, the total width of D S j (2317, 2460) is less than 100 keV. Both resonances are 
extremely narrow if they are cs states. A precise measurement of their total widths may help distinguish theoretical 
models of their quark content. 

As an straightforward extension, we also calculate the pionic widths of B S j mesons in the (0 + , 1 + ) doublet. With 
masses estimated in Ref. Q , their single and double pion decay widths are listed in Table IIIII 
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